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SUMMARY 

The stability of a rotating flow in a triaxial ellipsoidal shell with an imposed temperature 
difference between inner and outer boundaries is studied numerically. We demonstrate 
that (i) a stable temperature field encourages the tidal instability, (ii) the tidal instabil- 
ity can grow on a convective flow, which confirms its relevance to geo- and astrophysical 
contexts, and (iii) its growth rate decreases when the intensity of convection increases. 
Simple scaling laws characterizing the evolution of the heat flux based on a competition 
between viscous and thermal boundary layers are derived analytically and verified nu- 
merically. Our results confirm that thermal and tidal effects have to be simultaneously 
taken into account when studying geophysical and astrophysical flows. 

Key words: tides, tidal/elliptical instability, planetary cores, finite element numerical 
simulations, convection, stratification, Nusselt number. 



1 INTRODUCTION 



The tidal or elliptical instability comes from a triadic parametric resonance betw een two inertial waves of a rotating fluid and 
an imposed elliptic deform ation (see for instance the review by Kerswef 20021 . and references ther ein). It is a generic phe- 



nomenon m 



Leweke and Williamson 



turbulence (e.g.lWidna ll et al.lll974l; Moore and Saffmanffl975 ) and vortex dynamics (e.g. 

199g| ; Le Dizes and Laportell2002l ; iMeunier et ah 2002). There, the elliptical deformation comes from interactions between ad 
jacent vortices. But this instability is also studie d in geophysical and astrophysical ti dally deformed bodies: for example, its 
presence has been sugge sted (i) in binary stars ( Rieutord 20031 ; Le Bars et al. 2010l ) and accretion disks ( Goodman 19931 : 



Ryu and Goodmanl 1994|). wher e they could participate in the en ergy and angular momentum exchanges b etween neighbour 



ing systems I Lubow et al. 19931 ). and (ii) in planetary cores (e.g. Aldridj?e_i?t_al] 



1997 



play an important role in the inductio n of a magne tic field I KerswelHmdTlalkusf 199a : 
2009) and even in planetary dynamos (jrVlalkusll 19931 ) 



Cebron et al 



201 oh. where they could 



Lacaze et al. 



I I2OO6I : 



Herreman et al 



Previous studies of tidal instability have been perfor med for an isothermal fluid , starting from a solid body rotation 
or from an homogeneous isotropic turbulent rotating flow ( Fabiionas and Holm 2003j) . Nevertheless, in all natural systems, 
temperature differences are also present, which lead either to stratification or to convection. The previous studies need then to 
be reinvestigated in order to quantify the impact of a temperature field on the elliptical instability. Conversely, from a thermal 
point of view, many studies have been pe rformed regard ing the convective flow of an incompressible homogeneous fluid in 
a rotating spherical shell (see for instance I Aurnoul 120071 . and references therein). However, since most astrophysical bodies 
are tidally deformed, these studies have to be reinvestigated to account for the potential presence of an ell iptical instability. 



Note t hat a link between tidal resonance of gravito- inertial waves and thermal effects has been proposed bv lKumazawa et al 
(2006) n a geophysical context. However, this study considers the resonant forcing of a given wave by tides, which is different 



from the elliptical instability possibly induced by parametric resonance. Then, to the best of our know l edge, the only previous 
work which combines thermal and elliptical effects is the theoretical study of Le Bars and Le Dizesl (|2006h . who considered 
the influence of an established diffusive temperature profile. They then interpreted the thermo-elliptical instability in terms 
of gravito-inertial waves resonance and determined its linear growth rate by a local approach. In order to complete this work, 
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a guideline of our study is given by the following questions: (i) how is the growth rate of the tidal instability modified by an 
imposed temperature difference in the spheroidal geometry? (ii) Does the tidal instability grow over an established convective 
flow? (iii) Is the heat flux modified by the instability, and what are the scaling laws involved? In this paper, we tackle these 
questions using a systematic numerical study of the thermo-elliptical rotating flow in a fixed triaxial ellipsoidal shell. The 
numerical model is presented in section 2. The variations of the heat flux and the growth rate are systematically quantified 
and interpreted in section 3. Finally, orders of magnitude and conclusions are drawn for geo- and astrophysical flows in section 
4. 



2 NUMERICAL MODEL AND ITS VALIDATION 

As sketched in figure[TJ we consider the rotating flow inside a fixed triaxial ellipsoid of axes (a, b, c), with a constant tangential 
velocity U \J\ — {z/cj 2 imposed all along the outer boundary in each plane of coordinate z perpendicular to the rotation axis 
{Oz), where U is the imposed boundary velocity at the equator. An homothetic ellipsoidal core, in a ratio rj = 0.3, is placed at 
the center, with also an imposed constant tangential velocity r)U \J 1 — (z/c) 2 along its boundary. In the following, we use the 
mean equatorial radius R eq — (a + b)/2 as a lengthscale and we introduce the timescale by writing the tangential velocity 
along the deformed outer boundary at the equator U = Q,R eq . Then, the problem is fully described by the adimensionalized 
Navier-Stokes equations with four dimensionless numbers: the eccentricity e = (b 2 — a 2 ) /{a 2 + b 2 ), the inner to outer core 
ratio rj, the aspect ratio c/b, and the Reynolds number Re — Q R 2 q /u, or its inverse the Ekman number Ek, where u is 
the kinematic viscosity of the fluid. In all the computations presented here, we choose c equal to the mean equatorial radius 
Reg = (a + fo)/2, as in classical experimental setups (see e.g. Lacaze et al. ( 20051 ) ). This allows us to focus on the well-known 



stationary spin-over mode of the tidal instability. Note that for lots of tidally deformed planetary and stellar bodies, the c-axis 
is the shortest axis because of the rotational flattening, whereas the longest axis is the one pointing to wards the tide-raising 



body. I n this configuration, the first selected mode of the tidal instability close to threshold may change ([Kerswell and Malku; 



1 199Sl l. Cebron et al.l (^Plpl ) 1 ). giving rise to oscillatory motions. Nevertheless, we expect the main physical processes described 



here to remain valid for these other modes. 

In addition to these purely hydrodynamic considerations, we also solve the temperature equation with constant and 
uniform temperatures T and T e applied respectively at the inner and outer boundaries. To model the buoyancy, the density 
is assumed to vary linearly with the temperature and the Boussinesq approximation is used. Then, the gravity term of the 
Navier-Stokes equation is written: p g* = p(T e ) [1 — a(T» — T e )] g, , where p is the fluid density, a the coefficient of thermal 
volumic expansion, g* the local gravity which depends on the position and T* the local temperature. In our model, the 
gravity is calculated by solving the Poisson equation for the gravitational potential <f>, assuming that the outer boundary 
is an isopotential. Thus, the thermal part of the problem is controlled by two dimensionless numbers: the Prandtl number 
Pr = v I 're, where re is the thermal diffusivity, and the Rayleigh number Ra — (a go (Ti — T e ) d 3 )/(re v), where d — R eq {l — rf) is 
the shell thickness and go is the surface gravity at the outer boundary of the corresponding sphere, i.e. go — 4/3 n G R eq p(T e ) 
with G the gravitational constant. Note that the kinematic viscosity of the fluid is assumed independent of the temperature 
in this work. 

Finally, the problem solved is described by the following system of dimensionless equations: 

du 1 Ra 

- + u- Vu = -V P + - An - pr Re2 _ 



+ u . 7u = - 7fl+ _A u - ___ rgi (i) 



V ■ u = 0, (2) 

A <t> = 3, (4) 

g = -V<P, (5) 

where T is the dimensionless temperature anomaly T — (T» — T e )/(Ti — T e ). All boundary conditions have been described 
above. This numerical setup is a simple model for a tidally deformed liquid planetary core with a solid inner core and an 
imposed temperature contrast. Note that in a geophysical context, because of the compressibility, the temperature difference 
AT — Ti — Te to consider is the superadiabatic temperature difference between the inner core boundary and the core-mantle 
boundary (CMB). All dimensionless numbers as well as explored numerical ranges and typical values for the Earth's and Io's 
cores are given in table [T] 

Since there is no simple symmetry in the studied configuration, the fast and precise spectral methods classically used 
in planetary cores studies are not relevant here. Our computations are performed with a commercial software using a finite 
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element method, which allows to correctly reproduce the geometry and to simply impose the boundary conditions. An 
unstructured mesh with tetrahedral elements was created. The mesh element type employed is the standard Lagrange element 



P1 — P2, which is linear for the pressure field but quadratic for the velocity field. The temporal solver is IDA IjHindmarsh et al 



2005). based on backward differencing formulas. At each time step the system is solved with the sparse direct linear solver 



PARDISC0. 

The numerical procedure is the same for all the computations presented in this paper: starting from a fluid at rest in a 
spherical shell of radius R eq with a diffusive temperature profile, the shell is set in rotation at a constant angular velocity 
Q and a first statistically stationary state is obtained by solving the Navier-Stokes and temperature equations; then, the 
spherical shell is instantaneously deformed to the chosen triaxial ellipsoidal shape and we follow the dynamics. Note that this 
instantaneous deformation allows to study the growth rate of the elliptical instability on a well-known conductive or convective 
base flow. But we have also performed computations with a gradual increase of the deformation, as well as computations 
starting from a deformed shell with a fluid at rest before imposing the rotation at the boundaries: the classical isothermal 
tidal instability then grows, which is only slightly modified by then imposing a temperature difference. In all cases, the final 
steady-state is the same. 

Two quantities are systematically studied: (i) the growth rate a of the tidal instability, which quantifies the feedback 
of the temperature field on the elliptical instability, and (ii) the Nusselt number Nu, which quantifies the feedback of the 
elliptical instability on the thermal properties of the system. Since the elliptical instability induces a strong three-dimensional 
destabilization of the initial flow, the growth rate a is determined as the time constant of the exponential growth of the mean 
value of the vertical velocity W = -k f v \w\ dr, where w stands for the dimensionless vertical velocity and V for the volume of 
the ellipsoidal shell. The Nusselt number is calculated by dividing the measured mean heat flux at the outer boundary (once 
a statistically stationary state is reached) by the theoretical value of the conductive heat flux in a sphere of radius R eq . It has 
systematically been checked that the difference between this analytically known flux and the numerically calculated flux for 
a purely conductive ellipsoid is below 2% for our all simulations. 

The purely hydrodynamic part of the numerical simulation has already been validated in ICebron et al. I (|2O10l ). A com- 
plementary validation of the resolution of thermal effects has been done by computing the heat flux in a rot a ting s pherical 
shell: in figure [5] (a), our results at Re = 344 are compared with the compilation of results given in I Aurnoul (|2007l ) for the 
modified Nusselt number defined by Nu* = p r B , where E = n ^ is the Ekman number based on the gap width, i.e. 
E = ((l-r?) 2 Re)' 1 . Note that our numerical results fill an area wit h few data and agr ee with the general trend. Besides, our 
results are compared in figure [2] (b) with the scaling law proposed by Christensen ( 20021 ) obtained from numerical simulations: 

Nu* = 0.077 {Ra* Q f /9 (6) 

where Roq = Ra Nu E 3 /Pr 2 is the Rayleigh number based on the heat flux instead of the temperature contrast and which 
does not take into account any viscous or thermal diffusive effect. An excellent quantitative agreement is found. 



3 INTERACTION BETWEEN THE ELLIPTICAL INSTABILITY AND THE THERMAL FIELD 



In their analytical WKB ana lysis of the linear stability of elliptical streamlines in the presence of a diffusive temperature field, 
Le Bars and Le Dizesl (|2006l ) found an analytical expression of the inviscid growth rate of the elliptical instability in the limit 
of small deformations as: 



9 - 3Ra 
16 - 4Ra 



where Ra 



7)-l Ra E z 



(7) 

h — — — is a modified Rayleigh number. The expression reflects the modification by thermal effects of the 
elliptical inst ability, which in this c o ntext corresponds to the parametric resonance of 2 gravito-inertial waves with the tidal 
deformation. Le Bars and Le Dizes ( 20061 ) also showed that the elliptical instability cannot grow for Ra > 3 because no 
resonance between inertial modes is then possible. As shown in figure [3] (a), the same decreasing trend with Ra is recovered 
in the spherical shell, even if the results differ quantitatively from this theoretical local estimate. The first noticeable result 
from our complete calculations is that the gro wth rate of the elliptical inst ability is significantly enhanced by the thermal 
stratification (i.e. Ra < 0), as first suggested bv lLe Bars and Le Dizesl l|2006l ). It may seem at first surprising that an a priori 
stabilizing effect leads to an encouragment of flow destabilization, but such a behavior has already been shown for instance in 
a Taylor-Couette flow, where a stable configuration can be destabilized by the pr esence of an axial stratifica tion, leading to the 
so-called strato-rotational instability (see for instance the recent experiments bv lLe Bars and Le Gal 2007). The second result 
of planetary relevance is that the tidal instability can still grow on an established convective flow. As shown in figure |4j the 
organization of the flow is then completely changed. Indeed, in the absence of tidal instability, the convective flow organizes 
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in a classical columnar structure shown in figure [4] (a) (the so-called Busse columns), where motions are mainly bidimensional 
and vertical velocities remain small, coming only from the Ekman pumping inside the columns. On the contrary in figure[4](b), 
Busse columns are completely destroyed by the tidal instability, whose velocity field is fully three-dimensional. The growth rate 
a of the elliptical instability progressively decreases when the Rayleigh number increases as shown in figure [3] (a). Indeed, the 
thermal convective motions can be seen as a perturbation superimposed on the classical isothermal base flow of the elliptical 
instability, corresponding to a rotation along elliptical streamlines: when the Rayleigh number increases, the typical convective 
velocity increases and the ellipticity of the base flow streamlines is less and less felt by the fluid particles. The dimensionless 
version of this model can be quantified in defining the typical Rossby number of the flow, i.e. the ratio between the root 



mean square convective velocity and the typical rotating velocity £lR eq : Ro = y ^ J u • u dr, where u is the dimensionless 
velocity in the rotating frame. One would expect to find a critical Rossby number of order 1 for the disappearance of the tidal 
instability. As shown in figure [3] (b), the critical value Ro « 0.25, where the growth rate is zero, is indeed found for Pr — 1 
and Re = 344, which is near the threshold of the tidal instability. The critical Rossby number then seems to increase with the 
Reynolds number (see figure [3] (a) for Re = 688) , if it exists at all. Indeed, when the Rayleigh number increases, the typical 
convective velocity also increases whereas the typical lengthscale of convection L conv simultaneously decreases. In the limit of 
very large Rayleigh numbers, L conv is very small and turbulent motions become fully t ridimensional. We th e n exp ect thermal 
convective motions to be seen as an effective turbulent viscosity, where the results of IFabiionas and Holml l|2003h predict an 



encouragement of the tidal destabilization. The extensive study of this regime is very interesting, but currently out of reach 
of our numerical resources. Nevertheless, our current results allow to use the critical Rossby number Ro c = 0.25 determined 
previously for Re — 344 as a lower bound for flow at larger Reynolds numbers. This result is sufficient for the planetary 
applications that we are interested in, as will be shown in the next section. 

When thermal and tidal instabilities are simultaneously present, heat transfers are determined by the competition between 
the natural thermal convection and the forced convection due to the elliptical instability. In the latter case, the flow driven 
by the instability can be schematically described by three distinct regions: a fully mixed and isothermal bulk, and viscous 
boundary layers at the outer and inner boundaries, where the temperature is purely diffusive. The advective heat flux driven 
by the tidal instability is then determined by the size of the viscous boundary layers <5„, following Nu — 1 ~ d/8 v . With 
the classical scaling law of viscous boundary layers, we obtain Nu — 1 = ai/V~E*, where a?i is a constant, E* = - - d the 
Ekman num ber associated with the flow driven by the instability, and u e i the fluid velocity due to the elliptical instability. 
As shown in Cebron et al. (2010). this velocity u e ; may be computed numerically in the absence of temperature gradient by 



subtracting to the computed flow the theoretical base flow corresponding to an imposed rotation along elliptical streamlines. 
„"g generally scales as cti y/Re - Re c near the threshold of instability, where «2 ~ 0.036 has been determined numerically 
and where Re c is the critical Reynolds number for the onset of tidal instability equal to Re c = (5.24/e) 2 in the isothermal 
case for a stationary deformation. then saturates toward 1 far from the threshold. Keeping the same scaling in the 

presence of thermal effects, we expect the Nusselt number to vary as: 

Re N 1/4 



Nu - I = ai V<*2 (1 - n) - lj v^ite. (8) 

Our numerical results shown in figure[5](a) validate this scaling with «i w 0.01. Far from threshold, E* = E and the scaling 
law becomes 

Nu = ai/VE. (9) 

Note that according to this model, the Nusselt number is a measure of the amplitude of the instability. Then, the agreement 
between our results and the scaling law confirms that the amplitude of the instability is not inhibited by a strong stratification, 
provided that the Reynolds number is far enough from the threshold. 

The systematic evolution of the Nusselt number with the Rayleigh number is shown in figure[5](b). In the absence of tidal 
and convective instabilities, Nu remains equal to 1. In the presence of tidal instability, the heat flux remains constant at a value 
larger than 1 up to a transition Rayleigh number which depends on Re, where the natural convection becomes more efficient 
than the forced one. Indeed, (|8]) and (|9]) describe the heat flux by the forced convect ion only, and are valid in particular for 
stably-stratified regimes and before the onset of convection Ra < Ra c . But following King et al. 1 2009h . once Ra > Ra c and 



convective and tidal instabilities are competing, the heat transfer is determined by the most efficient mechanism. Following 
scaling laws (|6} and Q, the transition between the forced convection and the natural convection appears for the transition 
Nusselt number Nu* = N p E = 0.077 (Rag) 5 ^ 9 = an which corresponds to the transition Rayleigh number: 

Ra t ~2.5E- 8/5 Pr 1/5 . (10) 

For Ra > Rat, heat transfers are controlled by the natural convective motions, whereas for Ra < Rat, heat transfers are 
controlled by the tidal instability. For planetary applications, it is more convenient to write this result in term of flux Rayleigh 
number defined by Raf = Ra Nu, for which the transition appears at 

Ra ft « 0.025 E~ 21/1 ° Pr 1/5 . (11) 
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4 PLANETARY AND STELLAR APPLICATIONS 



As explained in details by ISumita and Yoshidal (120021 ) . a stable density stratification may have existed in the whole early 
Earth outer core and could even have been sufficiently large to prevent dynamo action. This initial stratification has been 
disrupted by some mech anism, which could have been either a gradual mechanism, including encroachment a nd entrainment 
I Lister and Buffett 199ct ). or catastrophic mechanis ms driven for instance b y the 1/3-annual forced nutation ( Williams! 1994 ) 
or by tidal resonant forcing of inertial-gravity waves ( Kumazawa et al. 20061 ) . Our work suggests an alternative proposal based 
on a parametric triadic resonance of gravito-inertial waves. Note that since the stratification enhances t he growth rate of th e 
tidal instability and since the early Earth was already unstable from a purely hydrodynamic analysis ( Cebron et al. 20101 1 ■ 
then a fortiori with a thermal stratification, the tidal instability can be a mechanism for this overturn, hence playing a role 
in one of the major events in geological Earth history. 

We can also apply the heat transfer scaling laws obtained in the previous section to the present convective state of the 



Earth's core. Following IChristensen and Aubert (|2006r ). the typical values are E « 5 



10" 



Ra, 



Q 



3- 10"" and Ra 



f » io 29 

0.41 



for Pr = 0.25. The first consequence th at we can deduce is that Ro ~ 10 according to the scaling law Ro — 0.85 (Rag) 
given bv lChristensen and Aubert |2006h , which is significantly below the lower bound for disappearance of elliptical instability 



Ro c « 0.25 found in section [3] That means that the convection does not prevent the development of the tidal instability. 
Besides, the scaling law (|11[1 gives the estimate Raj t ~ 2 x 10 28 . The Earth's core is just above the transition between the two 
kind of convections considered in section [3] and its heat flux is thus controlled by natural convection. However, the closeness 
of the values obtained here does not preclude a tidally dominated flux in previous times or in other planetary systems, and 
in any case highlights the relevance of the elliptical instability as one of the major process at work in planetary cores. 

To finish with, many authors have assumed that the top of t he Ea rth core is stably stratified, this layer being the so- 
called ocean of the core (see for instance IStanlev and Mohammadi (2008) for a recent review of the literature on this subject). 
Actually, as described in this recent work, many clues lead to the idea that several planets contain stably stratified layers in 
their electrically conducting regions. An opposite (but dynamically similar) situation occurs in stars under the photosphere, 
where the stratified radiative zone takes place under the so-called convective zone. Then, we can address the simple following 
questions: in such a bi-layer flow, can the elliptical instability grow? In which area, i.e. the stratified zone, the convective 
zone or both? And is the heat flux completely modified by the instability? A single simulation is presented in figure [6] as 
an illustration of the answer. In this simulation, the same temperature T = is imposed at boundaries, and a temperature 
T = 1 is imposed at an internal neutral boundary corresponding to an homothetic ellipsoid in a ratio 772 = 0.7 : this allows to 
obtain a stratified inner layer and a convective outer layer, with clearly visible convective cells (figure [6] (a)). Once the tidal 
instability is excited, figure [6] (b) shows that the instability grows over the whole fluid and controls the heat flux. Hence in 
presence of the elliptical instability, a thermal stratification (i.e. a subadiabatic temperature profile) is not synonymous of 
thermal isolation, which is a key point from a geophysical point of view. 



5 CONCLUSION 

This paper presents the first numerical study of the interaction between the tidal instability and a thermal field. Even if the 
range of parameters (e.g. the Reynolds number) accessible to our numerical tool is relatively limited, we have validated general 
physical processes of direct relevance for planetary dynamics. We have shown that a stable temperature field encourages the 
tidal instability. We have demonstrated also that the tidal instability can grow on a convective flow, and may disrupt the 
famous Busse columns in planetary cores. Finally, we have demonstrated that the heat flux at planetary scales may be 
controlled by the forced convection due to this tidal instability, which in any case plays a fundamental role in the organization 
of fluid motions. 
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Table 1. List of relevant dimcnsionlcss parameters with the explored numerical ranges as well as their typical values in the Earth's and 
Io's cores. 



Dimcnsionlcss parameters 

Eccentricity 

Core radius ratio 

Oblateness 

Reynolds number 

Ekman number 

Ekman number based on the gap 
Prandtl number 

Rayleigh number 
Flux-based Rayleigh number 
Nussclt number 
Modified Nussclt number 

Rossby number 



E 



Ra = 

Ra* Q 



Definition 

=■ - >> 2 -a 2 
V 

c/b 

Re = 

E k = 1 1 Re 
= ((l-r,) 2 Re)" 1 

Pr = u/k 

_ a gp (Tj-Te) d 3 

Ra Nu V E 3 /Pr 2 



n R 2 eq /u 



Nil Q total /Q conduction 



Nit* 



Ro -- 



Nu E 
Pr 



Numerical ranges Earth & Io's core estimates 



- 0.32 




10~ 7 - 10" 4 


- 0.3 




0.35 - ? 


0.86 - 1 




997/100 - 995/100 


< 10 3 




10 15 - 10 13 


> 0.001 




10" 1S - 10" 13 


> 0.002 




10" 1S - io- 13 


1 




0.25 - 0.25 


-3.8 • 10 5 - 


10 6 


10 23 - ? (< Ra c ) 


-2.6 • 10" 5 - 


0.24 


IO" 13 - ? 


1 - 5.6 




10 6 - ? 


> 0.002 




10~ 8 - ? 


- 0.4 




10~ 6 - ? 
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Figure 1. Schematic representation of the considered problem. The gravity field shown here in the equatorial plane is 
calculated explicitly, assuming that the outer boundary is an isopotential. 



Figure 2. Evolution of the modified Nusselt number in the case of a rotating spherical shell (e = 0) for Pr = 1, n = 0.3 
and Re = 344. (a) Our computations are compared with a compilation of nu merical and experim ental results given by Aurnou\ 



2007). (b) Our computations are compared with the scaling law 0) given b mChristensen S2002 ) 



Figure 3. (a) Influence of the Rayleigh number on the growth rate of the tidal instability for e = 0.317, c = P r = 1 
and Re = 344 or Re = 688. The growth rate given by equation ^ is also show n ; taki ng into account a classical damping by a 
surfacic viscous term —K/^Ek, where K is a constant of order 1 'Lacaze et al. 200d) . The values of K have been determined 
by matching the theoretical value with the numerical one at Ra = 0, giving K = 3.07 for Re = 344 and K = 2.91 for Re — 688. 
(b ) Evolution of the growth ra te of the tidal instability with the Rossby number, calculated following the scaling law given by 



Christensen and Auberi \200&) Ro = 0.85 {Roq) 0A1 . For Re = 344, the instability disappears for Ro ; 



0.25. 



Figure 4. Comparison between the convective flows in the absence of tidal instability (a, e — 0) and in the presence of 
tidal instability (b, e = 0.317,) for the same values of the Reynolds number Re — 344, the Rayleigh number Ra — 18762 and 
the Prandtl number Pr = 1. The vertical velocity is shown in three different planes. Also shown in (b) is an iso-surface of the 
norm of the velocity ||u|| = 0.19, which clearly exhibits the 'S' shape of the spin-over mode. 



Figure 5. (a) Influence of the Reynolds number on the Nusselt number for a convective flow with £ = 0.317, c = ^2^! 
Pr = 1 and Ra = 18762. (b) Evolution of the Nusselt number with the Rayleigh number in spherical and ellipsoidal (e = 0.317, 
c = ) shells for Pr — 1 a nd Re = 344 or R e = 68 8. Note that the classical dependance of the critical Rayleigh for the 
convection Ra c oc 7?e 4 ^ 3 (see Tilaner and Busst (1992) f or instance) can be observed: the decrease of the Reynods number 



modifies the threshold and then, increases the convective motions. 



Figure 6. Visualization of the reduced temperature T in the case of a stratified layer under a convective zone (Re = 344, 
c = ). The inner and outer boundaries verify T = and at a ratio 772 = 0.7, a temperature T = 1 is imposed such 
that the Rayleigh number based on the thickness of the shell is Ra = 18762, whereas the Rayleigh number based on 772 is 
equal to Ra2 = Ra (iL^i ~ 1512. (a) Temperature in the equatorial plane for a spherical shell of outer radius R eq — c. 
(b) Temperature in the equatorial plane and surface of iso-value |ju|| = 0.19 at saturation of the tidal instability for the 
corresponding ellipsoidal shell at e — 0.317. 
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Figure 1. 
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Figure 2. (a) Upper figure, (b) Lower figure. 
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Figure 3. (a) Upper figure, (b) Lower figure. 
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Max : 0.015 




Min : -0.25 

Figure 4. (a) Upper figure, (b) Lower figure. 
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Figure 5. (a) Upper figure, (b) Lower figure. 
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Figure 6. (a) Upper figure, (b) Lower figure. 



